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Abstract 



The construction of discrete scalar wave propagation equations in arbitrary inhomogeneous media was 
recently achieved by using elementary dynamical processes realizing a discrete counterpart of the Huygens 
principle. In this paper, we generalize this approach to spinor wave propagation. Although the construc- 
tion can be formulated on a discrete lattice of any dimension, for simplicity we focus on spinors living in 
1 + 1 space-time dimensions. The Dirac equation in the Majorana-Weyl representation is directly recovered 
by incorporating appropriate symmetries of the elementary processes. The Dirac equation in the standard 
representation is also obtained by using its relationship with the Majorana-Weyl representation. 

Resume 

La construction d'equations discretes pour des ondes scalaires se propageant dans un milieu arbitraire- 
ment heterogene a ete realisee recemment par I'introduction de processus dynamiques elementaires qui 
obeissent a un equivalent discret du principe de Huygens. Dans cet article, nous generalisons cette approche 
a la propagation d'ondes spinorielles. Bien que la construction puisse etre formulee sur un reseau discret de 
dimension quelconque, nous nous limitons au cas simple de spineurs definis sur un espace-temps de dimension 
1-1-1. L'equation de Dirac en representation de Majorana-Weyl est retrouvee en imposant aux processus 
elementaires certaines symetries appropriees. L'equation de Dirac en representation standard est egalement 
obtenue a partir de sa relation avec la representation de Majorana-Weyl. 

1 Introduction 

The formulation of wave propagation by using the Huygens principle was investigated some years ago by the 
Transmission Line Matrix Modeling method (TLM) This method retrieves the Maxwell equations by 

introducing current and voltage impulses which propagate along the bonds and are scattered on the nodes of a 
mesh of transmission lines. Refinements lead the TLM method to describe complex boundary conditions, gain 
or losses and propagation in inhomogeneous media for electromagnetic waves Such an idea was renewed 
recently for studying time-dependent wave propagation for scalar waves in inhomogeneous media |^ . Instead 
of considering currents or voltage pulses, the authors introduced some arbitrary scalar quantities that propagate 
on a Cartesian lattice. Both approaches include the action-by-proximity of the pulses, or of the scalar quantities, 
when they propagate from node to node in one time step and the emission of secondary waves at each node 
by means of a scattering process which emits the incident energy in all directions. Since voltage and current 
impulses are equivalent to electric and magnetic fields on a two-dimensional mesh, it is not surprising that the 
TLM method applies to the Maxwell equations. Nevertheless, the nature of the wave equations which result 
from such a formulation applied to scalar quantities propagating on a Cartesian lattice remains unclear. A 
beginning of answer arises in the work of Sornette et al. , who have enlarged the scalar model to construct 
the Klein-Gordon equation and the Schrodinger equation. However, the results were limited to homogeneous 
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media and the Schrodinger equation was ill-defined in the continuum limit. In order to generalize this attempt, 
a constructive approach was explored recently |^ to recover various kinds of scalar waves in an inhomogeneous 
medium. Starting from basic principles and incorporating fundamental symmetries to the scattering nodes, 
the authors derived in a systematic way time-dependent scalar wave equations in inhomogeneous media. They 
exhibited a unified equation which properly tuned by a unique parameter yields either the Klein-Gordon equation 
or the Schrodinger equation with a well defined continuum limit. This derivation offers a general framework, 
including the related TLM approach and opens up possible generalizations to describe spinor or vector wave 
propagation. 

In the present paper, we extend this previous work to the description of discrete wave propagation equations 
for free spinor fields. For simplicity, we restrict the presentation to the derivation of the Dirac equation on a one- 
dimensional regular lattice but the proof can be worked out on any underlying discrete lattice. Throughout the 
paper, we call currents the scalar quantities which obey a simple dynamics of propagation and scattering. The 
problem being linear, we define naturally the two scalar components of the spinor field as linear superpositions 
of those currents. In order to derive coupled propagation equations linking together both components of the 
spinor field, we first require closure conditions compelling the form of the matrices which describe the scattering 
processes. Finally, by taking into account appropriate space-time relativistic symmetries, the discretized Dirac 
equation in 1 -I- 1 dimensions is recovered. The current model is constructed for the two usual representations 
of the Dirac equation in 1 -I- 1 space-time dimension, i.e. the Majorana-Weyl representation and the standard 
representation. 

Independently from the above approach, different kinds of microscopic models describing spinor wave equa- 
tions have been proposed in the literature. Those models belong to two distinct classes : they are either based 
on a random walk description, leading to Euclidean-invariant spinor wave equations, or they rely on a complex 
hopping-type dynamics on a lattice, leading to relativistic-invariant spinor wave equations. The first class of 
microscopic probabilistic models deals with critical behavior of the free Majorana spinor field where the time 
variable t is turned into a parameter that controls the approach to a critical point (the system is at criticality 
when t oo). It is indeed well known |^ that local scalar field theories can be studied within a random 
walk representation, where the time t corresponds to the length of a brownian path. Among such approaches, 
Mc Keon et al. |9| have constructed a probabilistic model, based on binomial processes, to describe the spinor 
field in 1 -|- 1 space-time dimensions. As the Dirac equation is recovered from real stochastic processes, their 
method applies only to the real Majorana representation of the Dirac equation. Another random walk approach, 
including a spin factor, is developed in [^0|. The model is defined on three-dimensional lattices and exhibits 
different critical behavior depending on the value of the spin. The second class of models deals with the Dirac 
equation itself. The paradigm of those lattice models are the Susskind fermions [11| which have been extensively 
used in lattice gauge theories. A different but closed approach is exposed in [12|, where the model describes 
the dynamics of point-like spinless particles on a three-dimensional square lattice in the presence of complex 
hopping rates of modulus unity. Related works have been applied to study condensed matter topics such as the 
ground state of the Heisenberg antiferromagnet in two or three dimensions ]l3{ , or the energy spectrum of flux 
states JT^ ^Sj . The construction in terms of currents presented in this paper can be considered as belonging to 
that second class of models. 

The backbone of the paper are sections |^ and |^, where we describe the construction of the Dirac equation 



in the Majorana-Weyl and standard representations respectively. Sections 2.1 and 3.1 introduce the basis of 



the current model in both representations. Sections 2.2, 3.2, deal with the derivation of the discrete coupled 
propagation equations. The Majorana-Weyl construction turns out to be solvable by choosing suitable sym- 
metries in section 2^. The current model is achieved by computing the scattering matrices in section 2A. On 
the contrary, the standard representation cannot be constructed in the same direct way. Henc e, t o compute 
the scattering matrices we use the linear transformation linking both representations (section 3.3). Finally, 
section ^ is devoted to discuss further issues and applications of this work. 
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2 Majorana-Weyl representation 



2.1 Basic definitions : currents and fields 

This section is devoted to the construction of the Dirac equation in the Majorana-Weyl representation, which 
reads (see |l6| ) 



9* 



d_ 

dx 



mc 



-ai 



(1) 



where ^ — {ipL i'nY' the two-component spinor field and ct^, i = 1, 2, 3, denotes the Pauli matrices. 

The elementary bricks of the construction are called currents. The currents are complex numbers which 
propagate from node to node along the bonds of a Cartesian lattice, each bond carrying two currents propagating 
in opposite directions. For simplicity, the lattice is chosen to be one-dimensional and regular with a lattice mesh 
size a. At any time t, the system is completely defined by the values of all currents along the chain. The time 
variable is also discrete and r denotes the time unit. A current propagates between two neighboring nodes in 
one time step. All the currents propagate simultaneously, i.e. the outgoing currents, denoted Si, i = 1,2, leave 
the nodes at some time t and become incident currents on the neighboring nodes at time t + t. The incident 
currents are denoted Ei, i = 1,2 (Fig. [^). 
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Figure 1: Sketch of the propagation step. The outgoing currents Si, i — 1,2, on node n propagate in one time 
step to become incident currents Ej, j ~ 1, 2, on the neighbor nodes. 

After this propagation step, all the incident currents are instantaneously scattered. The scattering process 
is described by scattering matrices attached to each node. It transforms the incident currents Ei, i = 1,2, on 
one node into outgoing currents Sj, j = 1,2 on the same node (Fig. ||). 
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Figure 2: Sketch of the scattering step. The incident currents Ei, i = 1,2, on node n are scattered instanta- 
neously to become outgoing currents Sj, j=l,2. 

Since our aim is to describe the two chiral fields ij^L and V'iii two kinds of currents, L and R currents, 
are naturally introduced, leading to four currents propagating simultaneously on each bond. For clarity, the 
one-dimensional lattice can be pictured by two sublattices carrying the L and R currents respectively (Fig. ||). 
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Figure 3: The chain is spht into two subchains. The upper and lower subchains carry the L and R currents 
respectively. 

Due to the structure of the Dirac equation which displays a coupling term between the two chiral 
components ipL and ipn, we also introduce additional currents which propagate between the two sublattices in 
one time step as depicted in Figure ^. 
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Figure 4: Propagation step for the L and R currents. 

These additional currents, which are called the L or R node-current in the following, do not propagate along 
the chain but remain located on the same node. They can be considered as propagating along a internal time 
axis which couples the L and R sublattices at each node of the chain. Moreover, these node currents participate 
in the same scattering process as the propagating currents (Fig. H). 
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Figure 5: Scattering step for the L and R currents. 
Thus, the scattering process reads 



Si,L{n,t) = '^SijxEj^L{n,t), i= 1,2,3, 



(2) 
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S'fc,_R(n,t) = ^Ski,REi^]i{n,t), k = 1,2,3. (3) 
1=1 

where Sy,L and Ski.B. are the complex entries of the scattermg matrices Sl and Sr attached to each node. The 
ensemble of all the scattering matrices belonging to the chain defines the medium or the "background" in which 
the currents live. Once the scattering matrices are known, the evolution of the system at any time is completely 
determined given any initial values of the currents. Throughout the paper the scattering matrices are chosen 
to be identical from node to node indicating that the medium is uniform. However, the matrices Sl and Sr do 
not need to be identical. 

The dynamics includes all the features of a discrete counterpart of the well-known Huygens principle. Namely, 
the principle of action-by-proximity from one node to the nearest-neighboring nodes is taken into account by 
the propagation step, while the interferences and the radiation of the energy in all directions are taken into 
account by the scattering process. 

Let us introduce some convenient notations. The bonds attached to a given node are numbered as displayed 
in Figure ^ : 1 for the left bond, 2 for the right bond and 3 for the node bond. 

bond k bond k bond k bond k 



propagation ^ ' i 

* ll ^ node Hk 



U ml m ^ 



node n t) node node n ^^^(^n^t + t) Ej{nk,t + T) 

time t time t + t 

Figure 6: Notations for the bond linking the central node n to one of its neighbor node Uk and the propagating 
currents defined on this bond. 

Now consider a current Sk{n, t) leaving the node n along the bond k, k — 1,2,3 (Fig. |^). At the next time 
step, this current becomes an input current E on the neighbor node along the A:*'' direction. This neighbor node 
is labelled n^. Note that this notation implies — n. Moreover, viewed from node n^, the bond on which lies 
the input current E is denoted k, leading to labelling this current a,s E = E-i^{ni.,t + r). This notation implies 
1 = 2, 2=1,3 = 3. Hence, the relations linking the outgoing and the incident currents during the propagation 
step read 

Ek,L{n,t + T) = S^Jnk,t), fc = l,2, (4) 
Ek,R{n,t + T) = S-^j^{nk,t), fc=l,2, 

E3,Lin,t + T) = S3M{n,t), (5) 



E3,R{n,t + T) = S3^L{n,t). 



The propagation rules (||) take into account a transmutation of the node currents S'ax and S^^r into E^ r and 
Es.L respectively. Figure |^ summarizes the two steps defining the dynamics of the L, R currents. 

To pursue the construction of the Dirac equation within this two-species current model, the chiral fields need 
to be defined. The fields iJjl, iJjr are complex valued functions which are defined on the nodes of the chain. The 
linearity of the problem compels the choice for each chiral field V'l or -ijjR to be a linear superposition of the 
L or R currents respectively. As the outgoing currents can be expressed as a function of the incident currents 
according to the scattering equations (||) , (||) it is sufhcient to define the chiral fields in terms of the incident 
currents solely : 

3 

i>Liri,t) = ^Xk,L Ek,L{n,t), (6) 
fe=i 
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Figure 7: Scattering step (a) and propagation step (b) for the L and R currents. 



3 

ipR{n,t) = ^Xi^REi^ii{n,t), (7) 
1=1 

where Xk.L: a-re complex numbers. 

2.2 Discrete propagation equations linking ifj^ and i/jr 

Under certain assumptions, which we caU thereafter closure conditions, we show that the two chiral fields, ^l, 
ipR, satisfy the following discrete propagation equations 



'ipL{n,t + T) = fL[i^L{n\t'),i!R{n\t'),ipL{n" ,t" ,i!R(n" ,t"),. . .) , (8) 
^R{n,t + T) = fR{i^R{n',t'),ijL{n',t'),ijR{n",t"),Mn",t"),...), (9) 

where the fields, ipL, V'-R, on node n and at time t + r are functions of both chiral fields on the same node and/or 
on nodes n' , n", . . ., at previous times t', t", ... Equations (||), (^) are closed equations in the fields. This means 
that fL,R. are functions involving the two chiral fields only and do not depend explicitly on the incident currents 
of both kinds. Moreover, the number of terms in the right hand side of equations (^), (H) must be finite. Those 
two closure conditions strongly constrain the form of the scattering matrices Sl,r, so that one finds 



Sl = Pl +diag(^i,L,/X2,L,^3,L), (10) 
Sr = -Pr + diag(/ii,fl,^2,fl,M3,i?), (H) 



where 



{P)ij,a ^ Pi.aXj^a, = 1,2,3, a^R,L, (12) 

and Xj^a denotes the complex numbers appearing in the definitions of the chiral fields. The pi^a and pi^a are 
complex numbers to be determined. As usual diag denotes a diagonal matrix. The discussion leading to the 
form (|l0[), ( [ll| ) of the scattering matrices is exactly the same as for scalar fields Therefore, we have only 
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quoted the result here. According to equations ([lO| ) and (|Tl|), it is easy to check that the outgoing current Sk 
for each species sphts into a term proportional to the incident current Ek of the same kind on the same bond 
k and into a term proportional to the field of the same kind 



Sk,L = Pk,LlpL + fJ-k,LEk,L, 

Sk,R = Pk,RiJB. + fJ-k,REk,R, k = 1,2,3. 



(13) 
(14) 



As the constructions of the equations (^ obeyed by ipL and ipR are carried out in the same way, we only 
consider the equation for in the following. An analogous equation can be derived for the field ipn. We fix 
the space-time units, a and r, to unity. Let us consider the L field on node n at time t + 1 which is defined as 
(Eq. (i) 



i;L{n,t+l) = y" Xk,LEk^L(.n,t + l). 



k=l 



According to the propagation step, it can be expressed in terms of the outgoing currents at the previous time t 



ipLin,t + 1) 



fc=l,2 



^k,L S-^ ]^{nk,t) + A3,L S3^R{n, t). 



(15) 



In the remainder of this derivation, we will use drawings as often as possible since those representations are 
far more transparent than formal equations. For example, the definition of the L field (Eq. (H)), ij)L{n,t) — 
Yl\=i ^k,LEk,L{n,t), is represented schematically as sketched below 
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whereas equation (lla) is sketched 
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At time t, the outgoing bond currents appearing in equation ( |l5| ) were instantaneously scattered by Sl 
matrices. Hence, Si^l and S2,l are linear superpositions of L incident currents on nodes n + 1 and n — 1 
respectively, whereas 5*3, _r is scattered by an Sr matrix on node n. The resulting linear superpositions are given 
by the specific form (|l3|). Schematically, this scattering process is represented as follows 



2,L 



As./", S: 



3,L ^3,_R 



Efc=l,2 h,LPk L i'Lirik, t) + \3,LP3,R ^h{n, t) 



+ 



Ai,l(/«2,L -E'2,l) 



which reads analytically 



+ X! ^k,LPk^L^k,Li'^k,t) + X3,LP3..RE3,B.{n,t). (16) 
fc=l,2 

^ ^ ^ 

current term 

According to the propagation step, the current term of equation (^6|), which is the sum involving only the 
incident currents and not the chiral fields, is deduced from three outgoing currents at time t — 1. Those 
outgoing currents are all defined on node n, as sketched in Figure || 
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Figure 8: The upper drawing is a sketch of the current term in equation at time t, while the lower drawing 
is a sketch of the same currents at time t — 1 before the propagation step. 



which translates analytically through the following equation 

current term = ^ ^k,L^^-j:^L Skx{n,t - 1) + X3,LfJ.3,R S3,L{n,t - 1). 



(17) 



fc=l,2 



Again, the outgoing currents appearing in the right hand side of equation ( p^ are deduced from incident 
currents on the same node according to the scattering process (|3|). This leads to the following equation 



current term = ^ >'k,LPk,LfJ'k.L + '^^xPsx^M \ V'L(",i-l) 

+ X! ^k,LPk,LPk,L Ek,L(n,t - 1) + X3XP3,B.P'3.L E3^L{n,t - 1) . 



k=1.2 



new current term 

The new current term in the last equation is schematically depicted in Figure 
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Figure 9: Sketch of the new current term in equation (p^). 
Plugging the right hand side of equation (nsh into equation (Hq) leads to 



^pL{n,t + l) - ^ >^kXPkXPk^L + ^3,lP3,lP3.r\ i^L{n,t - I) 



, fc=l,2 
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k=l,2 

+ ^ {^^k,L^^k,L)>^k,L Ek,L{n,t - 1) + {n3^Rfl3^L)h.LE3^L{n,t - 1). (19) 
k=l,2 

Equation (^9|) is almost in the closed form we are looking for but still contains a current term. In order to get 
an equation which involves solely L, R fields without any L, R currents, it is sufficient to impose the following 
relations 

^J'k,L^J■kL = k = l,2, 

(20) 



where fj-L is a constant parameter. The relations ( pOD have been chosen in such a way that the last term of 
equation ( p^ ) becomes proportional to ipL{n, t — 1). The insertion of equations (|2^) into ( p^ ) leads to a discrete 
propagation equation of the type of equation (^ : 

Mn,t + l)-^iL^ 1^ Afc,LPfc,L \ ^ Afc,Lp^^VL(nfc,t) + A3,LP3,flV'fl,(n,t). (21) 

V fc=l ^/^'-f- / fe=i^2 

The discrete propagation equation obeyed by the R field is directly deduced from the equation for the L field 
by exchanging all the subscripts L and R : 

1pB.in, t + 1)- ^IR ( 1 + ^ ^k.RPkM \ ^^(-^^ t-l)= ^ >^k,R.Pk R i^Ri^k, t) + \3,RP3,L '^L{n, t), (22) 

V k=l '^'^^■R / fc=l,2 



where equation (|22| ) has been supplemented by sufficient conditions of the type of equations ( poj) 

Pk,RPk n = Mi?^, A: =1,2, 

(23) 

P3,LP3,R = Aifl • 

The comparison of the second equations in both systems (pO|), ( p3| ) leads to the equality : fj.^"^ = pr^ . In the 
rest of this section, we will denote this constant by /x^. Now, if the matrix parameters appearing in equations 
(pT]), (22) are chosen according to the following relations 



3 \ / 3 ^ 

^k,LPk,L \ 2 / 1 , ^k,RPk,R 



then a first-order discretized time derivative appears 



ipL{.n,t + l)-ipL{n,t-l) = ^ >^k,LPk,Li^Link,t) + X3^LP3,Ri^R{n,t), (25) 

fc=l,2 

■i/'_R(n,t + 1) - V'_R(n,t - 1) = ^ XkMPi:,^R'ipRi'rik,t) + \3,RP3.L'ipL{n,t). (26) 



fc=l,2 



On the same ground, the sums J2k=i 2 ^k,L,Rpj: RipL.Rink, t) entering equations (^5|), ( pq ) can be turned into 
first-order discretized spatial derivatives by an appropriate choice of the products Xk,L,RP^ r, for k = 1,2. 



With this assumption, the two coupled discretized wave propagation equations (25), ( pq ) can be identified with 
the Dirac equation in the Major ana- Weyl representation. However, we don't need this last assumption at this 
step of our derivation. In the next section we show that implementing some suitable symmetries on the scattering 
process constrains even more the L, R matrix coefficients and turns the terms J2k=i 2 ^k.L.RPk l R^^L.Rink, t) 
into first-order discretized spatial derivatives. 
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2.3 Symmetries of the scattering process 

Our goal is now to prove that a few well chosen assumptions on the symmetries of the scattering process leads 
to the unique determination of the L, R scattering matrices and to the discrete propagation equations (p5|), (|2^ ) 
in a form very close to the Dirac equation (|l|). The symmetries are implemented both on the currents and on 
the fields. In this model, the scattering process is described by matrices which are unitary due to time-reversal 
invariance and reciprocity. This property insures the local conservation, at each time step, of the "flux" of the 
L, R currents : (| 5*1 1 + | 52 | + | •S's | )l,r = (| -Ei | + | | + | £^3 I ) l,r, which is a stability condition on the 
dynamics of the L, R fields at large times. Finally, the last symmetry required to complete the construction is 
reflection symmetry [x t— > —x) which is known to connect the two chiral fields. 

2.3.1 Time-reversal invariance 

Time-reversal invariance for the currents can be formulated by stating that the two scattering processes, sketched 
in Figure |l^, are equivalent under a time-reversal transformation 




Figure 10: Time-reversal invariance of the scattering process. The scattering process for the L, R currents 
depicted in (a) is equivalent to the scattering process deduced from time- reversal symmetry (b). 

Note that the L and R currents are exchanged in this transformation. This corresponds to the exchange 
of the spinor components V'l and V'fi when considering the time-reversal transformation in the Dirac equation. 
This leads to the following relation between the scattering matrices 

Sl = {SR*r\ (27) 

According to the form (p^, (p]), ( |l^ of the L,R scattering matrices, equation leads to the following 
constraints on the L, R matrix coefficients 

C, 

for fc=l,2,3, (28) 

1 

C*' 

where C, C are two constants which obey the condition C + C* — ^k.LPl r- Equations ( ^^ allow to 

express the R matrix coefficients in terms of those of the L matrix. 



Pk,R 
Pk,L 

■^k.R 
^k,L 



C" 



fJ-k.R- 



^k,R 



^J'kM- 



Pk,L 
^ Pk,R 
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We must also consider time- reversal invariance for the chiral fields ipL and ipR- For its implementation, it is 
convenient to introduce the scalar field constructed as a linear superposition of outgoing currents 



ip^{n, t) = ^ KkSk{n, t), 



(29) 



fe=i 



where the coefficients are complex numbers. Up to now the scalar fields tp^, considered so far were linear 
superpositions of incident currents (equations (H), (^). Owing to the linearity of the scattering process (||), (^) 
there always exist numbers Kk such that the fields obey the same propagation equations as those satisfied 
by 1 namely equations (p5|), (p6|). By taking into account the special form ( p^ ) of the outgoing currents, 
equation (p9) for the L field becomes 



tpl{n,t) = '^Kk,LPk,L i>f{n,t) + y^^Kk,LfJ'k,LEk,L{n,t). 



\k=l 



k=l 



It is easy to show that the condition Kk,LlJ-k,L — i^L^k,Li where kl doesn't depend on the subscript k and is 
given by = 1 + X]fc=i ^k,LPk,L/ fJ-k,L, leads to the equality between the outgoing field and the incident field : 
ipf = tpf. Finally, using equations (M), V'f reads 



■4>f{n,t) = y^Kk,LSk,L{n,t) ^ fi^y^ JhLSk L{n,t). 

^ ^ tJ-k,L 



k=l 



k=l 



Similarly, one finds 



V'|(n, t) = ^ Hk,RSk,R{n, t) = p"^^ -^Sk,R{n, t) 



k=l 



k=l 



P-kM 



(30) 



(31) 



It is natural to implement time-reversal invariance on the chiral fields by the following conditions 



V^f(n,i) 
V'|(n,t) 



^|*(n,i), 



(32) 
(33) 



Here, we have explicitly distinguished the fields constructed on the outgoing currents, V'f rj and the fields 
constructed on the incident currents, i^f^ j^- Moreover the complex conjugate fields have been considered in the 
inverse process. Now if we write down the currents in equations (|3^), ( ^3| ) according to the definitions of the 
fields (H), (0), ( |29| ) the time-reversal invariance of the scattering process (Fig. |l^) leads to the following relations 



Kfe,L = K,R^ i^k,R = K,L for = 1, 2, 3. 
Then, using (pO|), ( [2^ ) with /i| = = /i^, the last equation (34) gives 

= 1, fc = l,2,3, and = 1, 

= /^^M^i?' k = l,2. 
Equations (^, (|34|) and (^ lead to : C" = /i^. Lastly, the relations (|8|) read 




P'kX- 



^kM 
Xk,L 

P*k,R 
Pk,L 



Pk,R 



IJ-k.R 







Xk,R 




PI,L 


1 


Pk,R 


c* 



for fc = l,2,3. 



(34) 



(35) 



(36) 



with C + = j:l^,Xk,L Pin- 
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2.3.2 Reciprocity principle 



The reciprocity principle states that the response of a system at a point r' due to an excitation at point r is 
identical to the reciprocal process, namely the response of the system at point r due to an excitation at point 
r' . In our current model, the application of this principle to the scattering process means that any scattering 
channel is identical to its reciprocal channel. The definition of a scattering channel is illustrated in Figure |ll| 
(a). An incident L current at node rt, i?3,L in Figure 111] (a), is scattered in three outgoing currents Si^h- 
A scattering channel is any of the three elementary processes, £^3,l — > E^^l S2,l and i?3,L — > S-i,L- 

Let us consider the first channel E^^l ^ S'i.l- The reciprocity principle states that this scattering channel is 
equivalent to the reciprocal channel for the R current illustrated in Figure |l^ (b). Roughly speaking, each 
scattering channel for L currents in one direction is equivalent to the same scattering channel for R currents in 
the opposite direction. As for the time-reversed processes considered in the previous section, note that the L 
and R currents are exchanged in the two reciprocal processes. This exchange of L and R currents is suggested 
by the exchange of the spinor components ■0L and tpn in the symmetries underlying the Dirac equation. 



Si^L 





Figure 11: Reciprocity principle of the scattering process. The scattering channel 3^1 for the L currents (a) 
is equivalent to the scattering channel 1 — > 3 for the R currents (b). 

Finally, reciprocity as defined above leads to the following relation between the scattering matrices 

Sl = '^Sr, (37) 
and to the following constraints on the L, R matrix coefficients 

Pk,L = l^kMi 

Pk,R = l>^kx, for fc = l,2,3, (38) 

A*fe,L = Pk,R, 

where 7 is a complex number. 

2.3.3 Unitarity of the scattering process 

The combination of equations (^^ and ( ^) shows that Sl and Sr are unitary matrices, each one satisfying : 

SS'^ = Id, (39) 
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where Id denotes the identity matrix. Unitarity which is expressed by combining together equations (pq), ( pq ) 
and (^ ) leads to 

\^^kx I = lMfe,-R, I = 1, 



A 



Pk,L 
Pk,R 



C_ 



for fc = 1,2,3, 



(40) 



where C/"f + (C/7)* = —A, with A = J2k=i 2 3 \-^k.L \ — J2k=i 2 3 l^fe,^? I ■ FinaUy, some manipulations using 
equations (p^, ( |3^ ) and yields the parametrization of the matrix coefficients in terms of Afe,L,Mfe,L 
and /i 



A. 



k,R 



Pk,R 



A \ ^ 
1 ( Pk.L 



^k,Lj 



A 



■^k.Li 



for A: = l,2,3, 



(41) 



Using this parametrization and choosing /i3^i ~ 1^3. r. = P from equations (|2C|), ( |3§| ) the scattering matrices iSl 
and S'fl, ((pi), (lu]), (P)) take the form 



SkLL 



Sij,R 



A \ ^ 



A 



^k,L'^l,L + P'k,LSkl, 



and the wave propagation equations (25), (pq) become 



(42) 
(43) 



tl;L{n,t+l)-^L{n,t-l) = 



A 



At,z,A2,LVL(n + l,i) + 



ipR{n, t+1) - -tJjR{n, t-l) = 



A 



M2,L 



Ai,LA^,iVfl(»^ + l,t) + 



A 



'iJjL{n,t). 



M2.L 



^i.LK,L'4'L{n- l,t) 



AT lA2,l Vfl(" - l,i) 



+ 7^3,^ i^Rin,t), 



(44) 



(45) 



2.3.4 Reflection invar iance 

Reflection invariance states that the Dirac equation (|l|) is invariant by the simultaneous exchange of V'l into 
tpn and X into —x. We transpose this property to the scattering process of the currents. Reflection invariance 
of the scattering process can be decomposed in three elementary steps sketched in Figures |12|, ^ and |l|. 
This symmetry is formulated by stating that the two scattering processes represented in figures (a) and (b) 
respectively, are equivalent under a reflection 

For example, the equivalence of the scattering processes (a) and (b) depicted in Figure n3 leads to the 



following relations between the matrix elements : Ski.L 



^kLR 



, fc = 1, 2, 3, and 1 = 1. Similar relations between 
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El, 



Si, 



2,L 



n — I 




n 


n + l 


n — I In 

S3,L* 


n + 


















scattering 


; 







(a) 








scattering 



ql 

'~3,R 
'~1,R - 





= S3,L 
S2,L 



S'i.R — Si, 



(b) 



Figure 12: Reflection invariance of the scattering process for the left bond current. The scattering process for 
the L current Ei l (a) is equivalent to the scattering process after a reflection symmetry for the R current E'2 
(b). 

Ski,L and s-^j ^ hold from Figures ^ and Finally, the constraints imposed by the reflection invariance on 
the scattering coefficients are summed up by the equation 

Sfci,L = Sfc,^^, for fc,; = l,2,3. (46) 
Taking into account the form of the L and R matrix elements ((p2[), (p3|)), we deduce from ( p6| ) 



|Ai,L 



= 1, 



(47) 



Ai,. 



A, 



3,L 



where, using ( p5D and the third equation of system (|38|), /ii.l and /i2,L can be written 



±1. 



(48) 



An additional constraint arises by requiring the reflection invariance of the propagation equations (44), ( |45| ) 
satisfled by the chiral fields. If we plug the last result (HS) in those two equations, the L, R fields are parity 



invariant provided the last two terms of the r.h.s. of equations ( [44| ) and (45) are equal, which implies 



7 = 



3,L 



e-y = ±1. 



(49) 



Using ( ^ ) and (|4g|), we obtain the phase of Xs^l 



'3,L 



(50) 
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2,L 



Si, 



2,L 



n — I 


n 


n+l 


n — I In 









S3,L * 













scattering 


i 














n + l 



(a) 








scattering 




(b) 



Figure 13: Reflection invariance of the scattering process for the right bond current. The scattering process for 
the L current E-) l (a) is equivalent to the scattering process after a reflection symmetry for the R current E'^ ^ 
(b). 



2.4 Identification with the Dirac equation and exact form of the scattering ma- 
trices 

In this section, we show that the propagation equations (jij), ( ^5| ) take either the form of the Dirac equation in 
the Majorana-Weyl representation or can be identified with two coupled Schrodinger equations (sec Appendix). 
We first note that the propagation equations (H), p5|) depend only on the ratios 



A 



2.L 



Ai, 



and V 



3.L 



A 



1,L 



(51) 



This is not surprising since the problem is linear which supposes that the chiral fields are defined up to a 
multiplicative constant : Ai^^. Furthermore, according to the equality between the modulus of A2,l smd Aix 
(equation (E^)) and to the definition of A, those two equations can be recast in the form 



A2,L 



ipji{n,t + 1) - ipBin,t 
(m* - / X2.L 



Ai, 



■ iJLin - l,t) 




\1,L 



Al,L 
A2,L 



^Rin,t), (52) 



(53) 



If we restrict the ratio (A2,l/Ai,l)*/(A2,l/Ai,l) to be equal to ±1, i.e. A2.L/A1.L is either a real number or 
a purely imaginary number, then the r.h.s. of the propagation equations (|5^), ( |53| ) displays a second order 
discretized spatial derivative or a first-order discretized spatial derivative, respectively. We postpone the case 
where A2,l/Ai,l is a real number to the appendix and concentrate on the other alternative where A2,l/Ai^i, 
is a purely imaginary number. Combining (A2.l/Ai_l)* = ^ (A2,l/Ai.l) with |A2,l/Ai.l| = 1 (Eq. (^7|)) leads 
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Si,L S: 



2,L 



n — I in 
Es,l\ 







n + l 



n — 1 



n + l 



scattering 



S3,L 





(a) 









scattering 



S3,R — ^3,11 



(b) 



^2.R — ^i, 



Figure 14: Reflection invariance of the scattering process for the node current. The scattering process for E^ l 
(a) is equivalent to the scattering process after a reflection symmetry for E'^ ^ (b) . 

to A2,l/Ai^l = i£x, where s\ ~ ±1. Writing [i — exp(i6') (see Eq. (p5[)), the propagation equations (|5^), ( |5^ ) 
depend only on the two signs e = eAEp, £7, and on the two real parameters \ v\ ^ 9. The space-time units, a and 
T, are again introduced so that cq — a/r is a. microscopic velocity. The choices e = — 1 and = 1 leads to the 
identification with the Dirac equation in the Majorana-Weyl representation (|l|) 

■^{ipL{n,t + T)-'iljL{n,t-T)) = 




a f ipLin + a,t) — ipL{n — a, i) ^ i j sin^j v |^ 

T2" 



2a J T \2+\p 



^Rin,t), (54) 



Zt 



2 sinO \ a f V'fi("- + a,t) — ipR{n — a, t) 




\2+\i^\yT\ 2a 
Comparison between (^), (^5|) and (|^) provides the velocity c and the mass m of the chiral fields 



Mn,t). (55) 



2sin6' , , 

Co— (56) 

1 sm0\uf 



(57) 



Exact form of the scattering matrices 



The L,R matrix coefficients (^2|), ( [43| ) depend only on /i, on the ratios /ill/m = M2,l/m = £/ii ^2,l/^i.l = 
ie\ = —iSfi, and i/. However the phase of v is fixed because v* /ly = i, thanks to the third equation of (4^ and 
to (|4^). Therefore, the L, R scattering matrices are parametrized by three real parameters : e^, 9 and 
introduce a complex number a defined by 



We 



a = ie^, (58) 
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so that the scattering matrices read 



St 



2 sin 6* 



Sr = - 



2 sin 6* 

2 + 



V 

Lastly, one finds easily the following relations 
and 



1 




V 




1 


2 


a v 


J/ 




1 


a2 


2 




1 


V 


V 


2 





Mdiag(£p,£p, 1) , 



(59) 



^diag(e^,e^, 1) 



(60) 



(61) 



(62) 



which will be used later in the standard representation (section 3.3) 



3 Standard representation 

In this section, the ideas developed in the construction of the current model for the Majorana-Weyl repre- 
sentation are reproduced for the standard representation in 1 + 1 space-time dimensions. Let us recall the 
Hamiltonian form of the Dirac equation in the standard representation 



9* 



d 

-T-cai — 
ox 



-0-3 



(63) 



where ^ — (V'+ ip-)'^ is the two component spinor field in the standard representation. Surprisingly, it 
turns out that the scattering matrices cannot be computed by a direct construction as in the Majorana-Weyl 
representation. More precisely, it appears that the natural symmetries of the Dirac equation in the standard 
representation transposed to the scattering process lead to contradictions in the computation of the scattering 
matrices. However, if we take advantage of its linear connection with the Majorana-Weyl representation, we 
are able to build indirectly the current model for the standard representation. 



3.1 Basic definitions : currents and fields 

The two current models constructed for the standard and the Majorana-Weyl representations share many 
features. Indeed, two kinds of incident and outgoing currents are defined on the bonds and on the nodes of a 
regular chain as the goal is to describe the propagation equations obeyed by the two scalar fields "0+ and tJj-. 
The lattice mesh size a and the elementary time step r are taken equal to unity. Transposing the notations of 
section 2.1, incident currents are denoted by -Efe,+, whereas outgoing currents are denoted by Sk'^+, Si'^-. 

Three incident currents of each kind live on each node of the chain depicted in Figure |l^. i?i,2,± denotes the ± 
bond-currents and /?3_± denotes the node-currents. The chain schematically represented in Figure |l^ is deduced 
from the chain depicted in Figure ^, for the Majorana-Weyl representation, by replacing the L currents by -|- 
currents and the R currents by — currents. 

The dynamics of the ± currents is similar to the dynamics defined for the L, R currents in the Majorana- 
Weyl representation. The ± currents obey a discrete Huygens principle on the chain which decomposes in two 
steps : a propagation step and a scattering step. At a given time t, all the ± currents undergo the dynamics 
simultaneously. The incident currents -E;.-, are scattered instantaneously on each node of the chain to 

become outgoing currents 5";,+ , Sk.-, which are linear superposition of those incident currents 



18 



l'--^/ n'':^/ + 1 



Figure 15: Sketch of the chain in the standard representation. The upper and lower subchains carry the + and 
— currents respectively. 



Si,+ {n,t) = ^s^j^+Ej^+(n,t), i = l,2,3, (64) 

3 

Sk,-{n,t) = J2''ki,-Ei,-{n,t), k = 1,2,3. (65) 
1=1 

where Sij,+ , Ski,^ are the complex elements of the scattering matrices S± attached to each node of the chain. 
Then, in a unit time step, the ± outgoing currents propagate from one node to the nearest-neighbor ones to 
become incident currents at time t + 1. Using the notations previously introduced in Figure ^, the propagation 
rules read 

Ek,+ {n,t + l) = Sj:_ink,t), fc = l,2, E3.+ {n,t + 1) = S3,+ {n,t), 

(66) 

EL^{n,t+l) = Sj^^{ni,t), 1 = 1,2, E3,^{n,t + 1) = -S3,-{n,t). 

Contrary to the construction of the current-model in the Majorana-Weyl representation (Fig. 0) the trans- 
mutation affects solely the bond-currents and not the node-currents. The choice of such propagation rules stems 
from the structure of the Dirac equation ( |63| ) in which the spatial derivative aid/dx couples and tp- while 
the last term (mc^ / fija^ does not. Additionally, a minus sign affects the propagation of the — node-current. 



The reason for introducing this minus sign will be clarified at the end of the construction (see section 3.2). The 
dynamics of the ± currents is summarized in Figure |l^. 

Finally, similarly to the definitions (^), (|^) of the chiral fields in the Majorana-Weyl representation, we 
suppose that the scalar fields ■0+ and ip- are defined on each site of the chain as complex linear superpositions 
of + or — incident currents respectively 



4'+{n,t) = ^Afc,+ i;fc,+ (n,i), (67) 

k=l 
3 

^-{n,t) = ^A,,_i;,,_(n,i), 



1=1 



where A^.^, A; _ are complex numbers. 
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\E2,+ / 



\E3 



time t 



S: 



3,+ 



scattering 



S3,- /■ 

time t 



(a) 



S: 



3,+ 



5* 



propagation 



2,+ 



3,+ 



(b) 



time t 



time t + 1 



Figure 16: Scattering step (a) and propagation step (b) for the + and — currents. 
3.2 Discrete propagation equations linking tp^ and tp_ 

Given the local rules which govern the dynamics of the ± currents, our aim is to derive two closed propagation 
equations involving the fields iIj± 



V'- (n, t + 1) = /_ (V.- in', t'), 4,+ {n', t'), (n", i"), ^+ {n" , t"), . . .) , 



(69) 
(70) 



where the fields V'i on node n and at time t+1 are functions of both fields on the same node and/or on 
neighboring nodes n', n", . . ., at previous times t', t", ... The necessary conditions constraining the scattering 
matrices S± for the equations (|69|), ( [70| ) to be closed remain the same as in the Major ana- Weyl construction 
(see section 2.2). So that the scattering matrices are of the form 



•5+ = -P+ +diag(/ii^+,/i2,+,Ai3,+), 
S- = -P- + diag(^i _,/i2 -,/i3 -), 



(71) 
(72) 



where P± are matrices which are generically given, for each kind of current, by (|l^ ) (with a — ±). The above 
form of the scattering matrices 5+ and S- implies for the outgoing currents the following equations 



Sk- = pk-ip-+Hk.~Ek-, k = 1,2,3. 



(73) 
(74) 
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The rest of the derivation is almost the same as for the propagation equations (^l]), ([2^) in the Majorana-Weyl 
representation. 



Using the propagation rules (66) and the definition (|67|), the field ^+ on node n at time t+1 can be expressed 



in terms of ± outgoing currents, defined at the previous time t and on the neighbor nodes Uk, 

^+(n,i + l)= ^ Afe,+ 5^_(nfc,t) + A3.+ 53.+ (n,i). (75) 



fc=l,2 

Equation (|75|) is sketched below 



n — 1. n '•. n + I 

Ai^_i_ S'2,- A2,+ S*!^- 



The outgoing currents entering equation ([75|) were scattered instantaneously at time t : Si^-, S2,-, and 53,+ , 
result from linear combinations of ingoing currents of the form (^3|) and ([74[). As a consequence, ip+{n^ t + 1) is 
represented schematically in terms of ± ingoing currents at time t as depicted as follows : 
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+ 



A3,+ (/^3,+ £"3,+ ) 



t 



which reads analytically 



k=l,2 



k=1.2 



current term 



(76) 



In virtue of the propagation and transmutation rules, the current term of equation ([76|) is deduced from three 
± outgoing currents that were defined on node n and at the previous time t — 1 



current term = ^ Afc^+^^ _ 5fc,+ (n, i - 1) + A3,+^3,+ 53,+ (n, t - 1). 



(77) 



k=l,2 



Applying the scattering rules to the outgoing currents entering the last equation gives rise to a field V'+ and to 
three + incident currents, at time t — 1 



current term = ^ Xk,+fJ.k,^Pk,+ + h,+fJ-3,+P3,+ \ ip+{n,t-l) 
+ ^ Xk,+H-i: _tik.+ Ek^+{n,t - 1) + X3,+fil _^_E3^+in,t - 1). 



(78) 



fc=1.2 



At this stage, we see that in order to close equation (|76|) at time i — 1 , it is sufficient to impose 



pi, A: = 1,2, 



= A* 



(79) 



where we have introduced the constant fj,'^. We end up the derivation of the closed propagation equation 
satisfied by ip+ by plugging the current term ( |7^ ) into equation (|76|), supplied with the two conditions (79), 
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X! -^fc.+Pfc - V'-(»^fc,i) + A3,+p3,+ V+(»^,i)- (80) 



fc=l,2 



The propagation equation for ■(/'_ can be directly deduced from equation ( pO[) by inverting all the subscripts + 
and — , except that we must take care of the minus sign involved in the propagation step ( |66|) of the — node 
current. The propagation equation for tp- is obtained in its final form by noting that the minus sign appears 
first in the expression of ip^{n,t + 1) involving ± outgoing currents (equation analogous to equation ([75)) 



ip-{n,t + I) = ^ Xk-Sj^_^{nk,t)-X3,^S3-{n,t). 

1 1 o 



A;=l 

Then, according to the scattering step, the last term of this equation gives rise to the term : — Aa^-pa^- ip- (n, t), 
which is analogous to the term A3. +/93. + 'i/;+(n, t) in equation (^0|). Finally, one can check that this term is the 
only one affected by the minus sign so that the propagation equation for ?/;_ reads 

\ k=i / 

X! ^fc -Pfe,+ V'+(»^fc,i) - A3 _p3 _ V'-(n,t)- (81) 

k=l,2 

Equation ( |8l| ) is supplied with the following sufhcient conditions 

(82) 

2 2 

Comparison between (^) and (|8^) shows that y?^ = which will be noted /i'^ in the following. 

3.3 Tranformation equations from the Majorana-Weyl to the standard represen- 
tation 

Implementing symmetries in the current model enables one to determine uniquely the form of the scattering 
matrices and of the propagation equations. Additionally, it provides a stability condition on the dynamics of the 
currents, which is fulfilled if the scattering matrices are unitary. Although this method was applied successfully 



in order to build the current model for the Majorana-Weyl construction (sections ^.3| and 2.4), it turns out 
that it fails for the standard representation. We thus proceed by using an indirect method in obtaining the 
± scattering matrices that relies on the linear transformation between the Majorana-Weyl and the standard 
representations. Contrary to the direct method, we show that there exist several possible expressions for the 
scattering matrices. The linear correspondence between the Majorana-Weyl and the standard representation 
reads 



I I I (83) 

According to the definitions of the fields ■0i and -0/? (Eqs. (|^), (0)), 0+ and -0- in ( ^3|) become 

3 3 

i'+ = ^KLEk,L + ^KREk,R, (84) 
fe=i 

3 3 

0- = 51 Vi^fc.i Vi?^fe,ii- (85) 



fe=i fe=i 
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If we remember that the Xk,L and the XkM are hnked together by the relations (|62|), then it is possible to express 
'!/'+ and 'ijj- as a function of the Xk,L or of the Xk^R uniquely. Let us first put the emphasis on the L parameters 
rather than on the R ones. For this purpose, let us write (|8^), ( p5|) in the following form 



V'+ = ^ X! [(^'^■■^ ~ 'i£ti>^k,R){Ek,L + i£i,Ek,R) + {Xk^L + i£p.Xk,R){Ek,L ~ iSt^Ek^n)] 

k=l,2 

+ \ [(A3,L + A3,ii)(S3.L + ^3,^) + (A3,L - A3.ij)(S3.L - ^3,/?)] , (86) 

V'- = ^ X! [(^'^■-^ + i^i^^k,R){Ek,L + i£i^Ek,R) + {Xk,L - i£fiXk,R)iEk,L - i£tj.Ek,R)] 

k=l.2 

+ \ [(A3,L - A3,fl,)(S3,L + S3.fl) + (A3,L + A3,flj(S3,L - ^3./?,)] • (87) 

If we plug (^) into the above equations, we see that terms are cancelled in each equation so that and t/j- 
become 

1p+ = Xi^L{Ei^L-i£^Ei^R) + X2x{E2,L+i£tiE2,R) + X3^L{E3^L + E3^R), 
1p- = Xi^LiEi,L+i£t,Ei^R) + X2x{E2,L-i£tiE2,R) + X3^LiE3,L- E^^b). 

Now, identifying Xk.+ and Afe._ with Xk.L, k = 1,2,3, imposes the following relations between the ± incident 
currents and the L, R ones 

Ei,± ~ Eix T iSfxEi^R, 

E2,± - E2x±^£^J.E2,R. (88) 

E3,± — E3X i E3^R. 

Assuming that the same relationships hold between the ± and L, R outgoing currents, one finds that by substi- 
tuting the propagation rules (Q), (||) for the L,R currents in (|8^), the propagation rules ( |66|) for the ± currents 
are recovered. The scattering matrices S± are also deduced from (|88|), by using the expressions of the scattering 
matrices in the Majorana-Weyl representation (|5^), ( |60| ) leading to 

S± = Sl. (89) 
If instead, we put the emphasis on the R currents rather than on the L ones, the ± fields are obtained through 



the equations (86), (|87|) where the roles of the L and R subscripts are exchanged. If we plug (p2) in those two 



equations, V'± becomes 



-0+ = >^l,R{El^R + ief,Eix) + h,R{E2,R-i£p,E2,L) + X3^R{E3^R + E3^L), 
V'- = —Xi^r{Ei^r — iSf^Eix) — X2,r{E2M + 'i£tiE2,L) — X3^r{E3^r — i?3,L)- 

Then, by choosing Afc_+ — —Xk,- — Xk,R, k — 1,2,3, we fix the relations between the ± incident currents and 
the L, R ones 

Ei,± = Ei_R±iSf^Ei^L, 

E2,± = E2M.Ti£t.E2,L, (90) 

E3,± = E3^R± Es^L- 

Again, if we assume that the same relationships hold between the outgoing currents of both representations, 
the propagation rules of the standard representation are recovered and the ± scattering matrices read 

S± = Sr. (91) 

Finally, we present a more symmetric transformation between the two representations where the L and R 
currents are treated on the same footing. The starting point is to write the ± fields (0), ( ^5|) in the following 
form 
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2.R) 



il>- = aXi^L^a* Ei^L + aEi^ji 
where a is defined by — ie^ (Eq. (58)). The identities : a*Xi^L = ctXi^R, aXi^L 



-^3,l(-E'3,l + E^^jf), 



a.*X2,L{aE2.L + ct*E2,n) + a^X^.^Lia*^ E^.l + o^Ez.r), 



(92) 
(93) 



aA2,i = Oi*X2,R, which foUow from (62), insure that neither the L incident currents nor the R incident currents 
are preferred in equations (|9^), (|93|). 
The choice 



fixes p+ and p_ 



a*Ai,L, 
aAiL, 



A2.+ 
A2,- 



aA2,L, 
a*A2,L, 



A3.+ 

A3,- 



Pl,+ = apl,L, P2,+ = Oi P2,L, P3,+ 

Pl,- = P2,- = aP2,L, P3,- 

and the relations between the ± incident currents and the L, R ones 



A 



3,L, 



^3,L, 



*2 



(94) 



(95) 



Ei,+ ai;i_L + Ei,^ = a* Eix + aEi^R, 

E2,+ = a*E2,L + aE2,R, E2,- = aE2,L + a* E2M, (96) 

£^3,+ = £'3,L + £'3,i?, £-3,- = E^^L + a'^Es^n. 

Assuming that the transformation rules ( p6| ) link also the outgoing currents of each representation, we recover 
the propagation rules of the standard representation and derive the ± scattering matrices 



/ 1 



2sin^ 



1 



1 1 



\ 



pdiag(£^,e^, 1) , 



(97) 



5_ = 



2 sin 6* 



/ 1 



/idiag(ep,£^,l). 



(98) 



Within this symmetric transformation both scattering matrices are different from Sl and Sr. Meanwhile, each 
one satisfies the symmetries : S± — '^S±, = S^., and are indeed imitary matrices. The first property. 



namely S± = 
property, S^^ 



'^S±, can be related to a reciprocity invariance of the scattering process, whether the second 
= S^, can be related to a time reversal invariance of the scattering process. But we stress that 
neither this reciprocity property nor this time reversal symmetry are symmetries obeyed by the components V'+ 
and ip- of the spinor field in the standard representation. Hence, contrary to the Majorana-Weyl representation, 
the symmetries of the current scattering process are different from the symmetries underlying the Dirac equation 
in the standard representation. This explains why we did not succeed in the direct construction of the standard 
representation. 

Finally, using (H), (H), equations (|^), (|l|) become 



■^{tp+{n,t + T) -tl;+{n,t-T)) = 

I 2 sin 9 \ a f ■(/'_ {n + a, t) — V-"- (n — a, t) 



2+kr / r 



2a 




■ip+{n,t), 



(99) 



2t 



{ij^{n,t + t) - ij_{n,t - t)) = 
I 2 smO \ a f ^J+{n + a, i) — ijj+{n — a, t) 



2+kr/ r 



2a 




(100) 
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where 



^2 i + 



k=l 



^k-Pk- 



k=l 



2 sin 6* 
2 sin 6* 



(101) 



2+ k 



|2 ' 



^2, + Pl,- — ^2,-Pl,+ 

^1, + P2,- = ^1,~P2,+ 

By using the relations (p6|), (|57|) one can check that equations (|99|), (100) are the discretized version of the Dirac 



2 sin 6*^1 



A3, + P3,+ 

, (|3) one 

equation in the standard representation. 

4 Conclusion and perspectives 

In this article we have recovered the Dirac equation in 1 + 1 space-time dimensions from a microscopic current 
model based on the Huygens principle and constructed from suitable symmetries, thus extending to spinor 
waves a previous work describing the propagation of scalar waves in an arbitrary inhomogencous medium 
This is not surprising since the current model combines all the basic features of wave propagation and since the 
symmetries are those underlying the Dirac equation. The originality of the method stems from the systematic 
construction based on very simple dynamical rules. Moreover the equation derived in this paper is independent 
of the type of discrete lattice and could be derived on any graph. The method developed in this paper is 
very close to the philosophy sustaining the construction of various cellular automata which are routinely used 
to solved numerically hydrodynamic or phase transition problems. Indeed, even if we cannot assign a precise 
physical meaning to the current variables, we choose the dynamics of those basic variables according to the 
fundamental laws of wave propagation. Hopefuly, the current model for the Dirac equation could be used as a 
useful tool to study numerically the time-dependent spinor propagation in complex systems with any kind of 
boundary conditions. 

A natural extension of this work would be to describe spinor field propagation in inhomogencous media. 
Indeed, a wide range of condensed matter problems are related to the random mass Dirac equation in 1 -I- 1 
space-time dimensions and could be studied by using such a model. Random mass Dirac spinors show up in the 
problem of a one-dimensional metal with a half- filled electron band and random backscattering |p7| , spin-ladder 
or quasi-one-dimcnsional spin Peierls systems In order to derive a Dirac equation with a mass and a 

velocity varying with the space coordinate, we can use the method developed in ^j. The description of scalar 
waves in an inhomogencous medium was achieved with the help of additional currents that trap a fraction of 
the "energy" on each node. Such additional currents can also be introduced for the Dirac equation. 

Another extension of this work would be to construct the Dirac equation in higher space dimensions. This 
should be useful for instance to study the plateau transitions in the integer quantum Hall effect where the Dirac 
equation shows up in 2 -I- 1 space-time dimensions ]l9| . Though the calculations are rather cumbersome, the 
Dirac equation in 3 + 1 space-time dimensions can also be constructed with our method. 

Finally, it would be interesting to see whether the current model described in this paper could be useful as 
a starting point for a path- integral formulation of the usual Dirac equation in 1 + 1 space-time dimensions | |20[ | . 
Despite the fact that the formulation of path integrals for spinor fields in terms of non-commutating variables 
is not new within the framework of quantum field theory, this description suffers from a lack of physical 
interpretation. As an example the field-theoretical description of a spinor field doesn't describe the way the 
fermion particle moves along a selected path as sketched by Feynman Some progresses in this direction 

have been done by extending the analogy between quantum mechanics and brownian motion to the Dirac 
equation j2^, ^ ^ ^ . But this point of view focusses on a probabilistic interpretation of the Dirac particle 
through Poisson processes. On the contrary our model deals with the Dirac equation itself and should give a 
direct interpretation of the paths followed by a spinor particle between two points. 
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Appendix: Construction of the Schrodinger equation with a spin 



This appendix deals with the construction of a Schrodinger equation for a quantum particle with a spin 1 /2 in 
1 + 1 space-time dimensions. This equation reads 



ih 



5* 

dt 



2m dx^ 



Id + vui 



(Al) 



where — {ipL i^R)^ is a two-component field. In section 2.4, we have obtained two discretized propagation 
equations for the chiral fields (Eqs. (Q), (^5|)) where the order of the discrete spatial derivative is controlled 
by the ratio {\2.L / I {^2.L I ^i.l) ■ To derive the Dirac equation, we have chosen \2,l/^i.l to be imaginary 
in order to obtain a first order spatial derivative. If we assign A2,l/Ai,l to be real then a second order spatial 
derivative arises. We write \2x/^i,l = £\ where e\ = ±1. The propagation equations (p4), (pq) depend now 



on the two signs e = £a£/^, £7, and on the two real parameters 



Then, we recognize the discretized version 



of the Schrodinger equation with a spin (Al) if we fix £ = —1 in the two wave propagation equations 



— {ipL{n, t + r) -ipLin^t- r)) = 




/ V-'L(n + a,t) + 4'L{n - a,t) - 2ipL{n, t) 



, , , 1 2sm9 



(A2) 



2t 



(V'_R(n, t + r) - iJjRin, t - t)) ^ 

I smO \ a? f ijjjiin + a,t) + ^j]i{'n ^ a,t) ^ 2tljji{n,t) 



|2 \ 



^ \ 2+\v\ 



tpL{n,t) 



1 / 2sin( 



r \2+k| 



(A3) 



provided that 



h 

2m 



sm ( 



a 

T 



2 h' 



(A4) 



V 

h 



i0\iy\ 



11 
r 



(A5) 



The L,R scattering matrices are parametrized by three real parameters : £fj,, and ji'l, as v is a purely 
imaginary number (using (B^) and (48)). Introducing again a (Eq. (58)), the matrices read 



/ 1 



2 sin 6* 



la 



\ 



/zdiag(£^,£^,l), 



(A6) 



- 9 -92 
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Sfl — ; 2^^^ 



( 1 ic? ic?v \ 

2 sin 6* 



/idiag(e^,e^,l) . (A7) 
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